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Heavy polar molecules can be used to measure the electric dipole moment of the electron, 
which is a sensitive probe of physics beyond the Standard Model. The value is determined 
by measuring the precession of the molecule's spin in a plane perpendicular to an applied 
electric field. The longer this precession evolves coherently, the higher the precision of the 
measurement. For molecules in a trap, this coherence time could be very long indeed. We 
evaluate the sensitivity of an experiment where neutral molecules are trapped electrically, and 
compare this to an equivalent measurement in a molecular beam. We consider the use of 
a Stark decelerator to load the trap from a supersonic source, and calculate the deceleration 
efficiency for YbF molecules in both strong-field seeking and weak-field seeking states. With 
a 1 s holding time in the trap, the statistical sensitivity could be ten times higher than it is 
in the beam experiment, and this could improve by a further factor of five if the trap can be 
loaded from a source of larger emittance. We study some effects due to field inhomogeneity 
in the trap and find that rotation of the electric field direction, leading to an inhomogeneous 
geometric phase shift, is the primary obstacle to a sensitive trap-based measurement. 



1 Introduction 

The permanent electric dipole moment (edm) of the electron, or other fundamental particle, 
can only be non-zero if both parity (P) and time-reversal (T) invariance are violated |Iil|2l. 
The weak interaction violates P, while T violation is equivalent to CP violation provided CPT 
invariance is accepted. The CP violation of the Standard Model generates edms that are far too 
small to detect . It also fails to account for the predominance of matter over antimatter in 
the universe 14). Many of the proposed extensions of the Standard Model contain new sources 
of CP violation and result in vastly larger edms |i5)|6J; indeed, some models predict values for 
the electron edm that are very close to the current experimental limit. A measurement of a 
non-zero edm would be firm evidence for new physics. 

Heavy, paramagnetic atoms and molecules offer high sensitivity to the electron edm. In 
an applied electric field, Ea, the change in energy of such an atom or molecule resulting from 
the electron edm, de, is deP{Ea)Eint- Here, Ei^t is a structure-dependent effective electric 
field whose magnitude scales as the cube of the nuclear charge |7 |, and P{Ea) is the degree of 
polarization of the atom or molecule in the applied field. The most sensitive measurement to 
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date, de = (6.9 ± 7.4) x 10~^®ecm, was made using a beam of Tl atoms HI, for which Eiat 
is large but P is small for realistic laboratory fields. Because polar molecules can have much 
larger values of P they offer even higher sensitivity |9, 10l|TTl. Using a continuous, thermal 
beam of YbF, we made the first measurement of the electron edm with a molecular system IIT2I . 
By upgrading this experiment to a cold, pulsed source, its statistical sensitivity has been greatly 
improved, and a new measurement is underway jT3l [T4l . Using the most recent advances in the 
production and manipulation of cold molecules, it seems likely that still higher precision can 
be obtained. The ability to trap polar molecules electrically 1 15 1 is particularly attractive for 
precision measurements such as the electron edm, because of the very long interaction times 
that are then available. 

The shot-noise limit on the statistical error in a molecule-based measurement of the edm, 
Sde, in the usual units of e cm, is given by 

Sd,^^ ^ ^. (1) 

e iPlEi^tWN 

Here, i^int is measured in V/cm, r is the amount of time each molecule spends inside the 
apparatus (the coherence time), and N is the total number of molecules that participate in 
the experiment. The inefficiencies of state preparation and readout, and any intervening state 
manipulations, are incorporated in our definition of N. In our current edm experiment, pulses 
of YbF molecules created at 4 K in a supersonic source travel at 600 m/s through a 60 cm long 
interaction region where the applied electric field is 20 kV/cm. There are 25 pulses per second, 
approximately 2500 molecules are detected in each pulse, and 50% of the running time is 
devoted to data-taking. For YbF, Ei^t = 26 GV/cm, and the polarization factor at 20 kV/cm is 
|P| = 0.7. Using these values in Eq. ([T]l, we estimate our shot-noise limited sensitivity to be 
7 X 10~^^ e cm/-y/day. As is often the case, the experiment operates slightly above the shot 
noise limit, magnetic field noise and source intensity fluctuations being the dominant sources 
of additional noise. At present, our statistical error is about 1.3 times the shot-noise limit. 

To increase the sensitivity, we have to increase the product r. There are, of course, 
many possible ways to do this. An increase in N might be obtained by improving the detec- 
tion efficiency or by increasing the time-averaged flux of molecules from the source. These 
would benefit both the beam experiment and the trap experiments considered here. Another 
strategy is to increase the fraction of molecules transmitted from source to detector, which 
might be achieved using molecular optics near the source to collimate the beam. For that to 
yield an improvement, the transverse phase space acceptance of the optics must exceed that 
of the machine without the optics. We will discuss beam optics below in the context of Stark 
decelerators, and we will see that the acceptance of our current machine is akeady comparable 
to that of realistic beam optics. An increase in t could be achieved without slowing or trapping 
the molecules, simply by making the machine longer. For a freely diverging molecular beam 
however, there is no advantage to increasing the interaction length L since the coherence time 
is proportional to L and the number of molecules reaching the detector scales as L^^, leaving 
t\/N constant. If optics are used near the source to collimate or focus the molecular beam, 
then an increase in L can be beneficial. 

In this paper, motivated by the prospect of an enormous increase in coherence time, we 
consider the possibility of measuring the electron edm using electrically trapped molecules. 
Given the difficulties involved in producing trappable molecules, we can expect the increase in 
r to be accompanied by a reduction in N . We evaluate the efficiency of Stark deceleration and 
so obtain an estimate for the statistical sensitivity that could be obtained. We then consider a 
few of the very many systematic effects that may scupper the measurement. We focus mainly 
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Figure 1 : Electric field dependence of (a) the Stark shifts of {N = 0— 6, A/ = 0) states of YbF, 
(b) the minimum number of deceleration stages needed to stop 340 m/s YbF molecules in a 
selection of these states, and (c) the polarization factor of each state. Each curve is labelled by 
the value of the rotational quantum number, N. 



on a measurement using YbF molecules, though our considerations are easily modified to other 
relevant molecular systems. 

2 Deceleration 

If the molecules are to be trapped, they must have a low velocity. We consider the prospects 
of bringing to rest the YbF molecules produced by our current supersonic source lfT6l . We 
take 340 m/s and 4K as the starting speed and temperature, these being typical when using 
room temperature xenon as the carrier gas. We have obtained speeds as low as 290 m/s by 
cooling the xenon, but only at the expense of molecular flux. It seems possible that a buffer 
gas source ifTTl [Tsl will provide slower YbF beams with higher intensity, but this has not yet 
been demonstrated. 

In designing a suitable decelerator, we can choose which rotational state of the molecule to 
use. Figure[TJa) shows the Stark shifts of the {N, M = 0) rotational states of YbF in fields up to 
250 kV/cm. Here, N and M are the quantum numbers of the rotational angular momentum and 
its projection onto the electric field axis. In the ground state the molecule has a negative Stark 
shift and can only be decelerated using the alternating gradient (AG) method |[T9l |20l |2T1| . 
The other states shown have positive Stark shifts at low field, and negative Stark shifts at 
high field. Molecules prepared in these excited rotational states could be slowed down with a 
"conventional" Stark decelerator for weak-field seekers |22|, which we will refer to as a WF 
decelerator. It would have to be operated at fields below the turning point of the state. Figure 
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[Tp3) shows the minimum number of deceleration stages needed to decelerate YbF molecules 
from 340m/s to zero. Deceleration in the ground state (Sec. 2.2 1 is obviously advantageous in 
terms of stages needed, but the difficulty of implementing the AG method is a disadvantage. 
Deceleration in the weak-field seeking states (Sec. 2. 1 1 is easier to implement but requires more 
stages. The best choice of state depends on the operating field. At 200kV/crr{2]the number of 
WF deceleration stages is minimized for N ~ 5. Lower rotational states are not as efficient 
because their turning points are reached at lower fields, as illustrated by the curve for TV = 4 in 
the figure. Higher rotational states are also less efficient because their Stark shifts are smaller 
at this field, as illustrated by the curve for iV = 6. 

We have looked at the rotational state dependence of deceleration, but what about the edm 
experiment itself? Its sensitivity depends on the rotational state and the applied electric field 
through the polarization factor. This is given in terms of the Stark shift, W, the electric field, 
E, and the dipole moment, ^,hyP — — (l//i) dW/dE. Figure [T|c) shows this polarization 
factor as a function of electric field for the various rotational states with M — 0. At any given 
field the ground state offers the largest polarization factor, which is one reason why our current 
experiment uses this state. An experiment that used a weak-field seeking state would have a 
smaller, though still useful, value of \P\. For example, the N ^ 5 state discussed above has 
IPI =0.36at£; = 130kV/cm. 



2.1 Deceleration in a rotationally excited state 

In this section we consider a practical decelerator for molecules in the {N, M) = (5, 0) state. 
We consider a design with a 4 mm square aperture and a periodicity of 2L = 24 mm, similar 
to the one used in Berlin to decelerate OH ll23l . Using voltages of ±40 kV, and a synchronous 
phase angle of 66°, 444 stages are required to decelerate from 340 m/s to zero. Each unit cell 
of the decelerator has two deceleration stages, one that focusses along x, and the other along 
y, so the overall length is 5.328 m. These are the parameters we use in the rest of this section. 

To estimate how many molecules we can expect to decelerate, we need to calculate the 
phase-space acceptance of the decelerator We start with the longitudinal acceptance. Since 
for our case the Stark shift is not at all linear (see Fig.[T]l, we will need to take a short diversion 
to extend the usual treatment of the longitudinal acceptance e.g. Il24l|25]| . The change in kinetic 
energy of the synchronous molecule in each deceleration stage is usually well approximated 
by A/\ = — A/fniax sin0Oi with AXmax being the maximum possible energy change per 
stage (a positive number), and (/iq the synchronous phase angle which goes from to 2n over 
the period 2L. For our case, this is a poor approximation, and we find it necessary to retain 
an extra term in the Fourier expansion of /S.K. We shall derive a general formula for the 
longitudinal acceptance, and then apply it to our case. We write the change in kinetic energy 
using the expansion 

AX = Aif^ax sin(n0o). (2) 

n odd 

Here, the a„ are coefficients that depend on the details of the decelerator and Stark shift of 
the molecule, and satisfy X^nodd'^" = 1- Next, we use the approximation of a constantly 
acting force to write down the equation of motion of a non-synchronous molecule in terms of 
its relative phase, 0, and its relative velocity v. 



= Aif„iax V (-l)("+^)/^a„ [sin(n0o + n4>) - sin(n(/)o) 



n odd 



(3) 



*We have operated decelerators at this field. Higher fields are very difficult to sustain. 
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where m is the mass. Integrating this equation of motion we obtain the contours of constant 
(relative) energy, 



V 
Vi 



71 odd 



U"+l)/2, 



■ coa{n(/)Q + ncf)) + 0sin(n0o) 



constant, (4) 



where Vi is the initial velocity and A'min is the minimum number of deceleration stages needed 
to bring the molecules to rest. The potential well represented by the second term in Eq. Q has 
a maximum at = tt — 20o, iu'espective of the values of the a„, for any sensible decelerator 
configuration. This is easily verified by inspecting the derivative of this potential with respect 
to (p. A molecule that reaches this maximum with zero relative velocity is moving on the 
separatrix between bound and unbound motion. The equation of the separatrix is thus found to 
be 
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n odd 



(-l)("+i)/2c 



■ COS (^n(f>o + ^n(j)j cos (^^ncf^ 

+{2(f>Q + (j) — n) sm{n(f>Q) 



-0. (5) 



For our case, the first two terms in this sum give a sufficiently accurate result. Using the 
electric fields obtained from a finite element model, and the Stark shift shown in Fig.[TJa), we 
can determine A A' numerically. Fitting this to the first two terms of Eq. with 03 — 1 — ai, 
we obtain oi = 0.825. The bold line in Fig.[2|a,c) shows the longitudinal acceptance area 
calculated using the first two terms of Eq. (j5]l, and this value of ai. Integrating the area inside 
this separatrix gives a longitudinal phase-space acceptance of = 23.4mmm/s. Note that 
the calculated acceptance area is smaller, 15.2mmm/s, when only the first term in Eq. Q is 
used (i.e. ai = 1). 

Next, we estimate the transverse acceptance, following the procedure outlined in e.g. 1(191 
I27I . The restoring force in the transverse direction is approximately linear in the transverse 
coordinate, ~ —kx{z)x, but with a spring constant that is a strongly-varying function 
of the longitudinal coordinate, with period 2L. The wavelength of the transverse oscillation 
is very much longer than 2L, except when the molecules have been decelerated to very low 
speeds. We calculate a mean spring constant, k^, by integrating kx{(t)) from (po — n to 4>o in 
the first switch state, and then from (pQ to (pQ + n in the second switch state, finally dividing by 
27r. In this way, we find the mean transverse angular oscillation frequency for our molecules 
of mass m to be uj^ = \/kx/m = 27r x 154Hz when (po = 66°. The corresponding phase- 
space acceptance in the x-dkection is — ttuIxTq, where rg = 2 mm is half the size of 
the decelerator's aperture. Thus, this approximate calculation gives us A^ = 12.2mmm/s. 
The bold line in Fig.[2|b,d) shows the transverse acceptance obtained using this approximate 
method. The total 6D acceptance is A ^ A^AyA^ = AlA^ = 3487 mm^ (m/s)^. 

The acceptance of a real Stark decelerator tends to be less than the idealized estimates 
given above. Coupling of the longitudinal and transverse motions can lead to unstable regions 
in phase-space |27|. Further loss occurs at very low speeds as the wavelength of transverse 
oscillations becomes comparable to the decelerator periodicity and the molecules are thrown 
out of the decelerator in regions where the focussing is weak |28 1. To calculate the true phase- 
space acceptance, we simulated the motion of molecules through the decelerator. Figure |2] 
shows the results obtained for 200 stages, where the final speed is 252 m/s, and for 444 stages 
where the final speed is 14 m/s, slow enough that one final stage - the trap - will bring the 
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Figure 2: Phase-space acceptance of a Stark decelerator for YbF in the (5, 0) state. Dots 
show the initial positions in phase space of those molecules transmitted and decelerated, as 
determined from a three dimensional numerical simulation. Bold lines show the acceptances 
calculated analytically as discussed in the text, (a, c) Longitudinal acceptance after 200, 444 
stages, (b, d) Transverse acceptance after 200, 444 stages. 



molecules to rest. To produce these plots, the trajectories of molecules chosen at random from 
a larger phase-space volume were simulated through the decelerator. The initial coordinates of 
successful trajectories, those that reached the end and were decelerated, are marked by the dots. 
The acceptance is the volume of the initial distribution multiplied by the successful fraction. 
After 200 stages, the region of phase space that contributes to the slowed beam is very similar to 
that predicted by the simple theory outlined above. In the longitudinal direction, part (a), some 
molecules are accepted from regions that are slightly outside the bold line. This is because 
off-axis molecules that travel closer to the rods of the decelerator experience a slightly deeper 
longitudinal potential than those on the axis. The region inside the bold line is not uniformly 
filled, the area closer to the origin having a higher density because it is less susceptible to 
coupling of the transverse and longitudinal motions. In the transverse direction, part (b), we 
see again that a few molecules outside the region given by the simple theory contribute to 
the acceptance because the transverse confinement is greater for some of the non-synchronous 
molecules. Once again, the central region has a higher density of accepted particles. The total 
6D phase-space acceptance for 200 stages is 1 168 mm'^ (m/s)-^, which is a factor of 3 smaller 
than the idealized prediction. An increase in the number of stages to 444, so that the molecules 
can be trapped, reduces the acceptance to 352 mm^^ (m/s)''. In the longitudinal direction, 
Fig.|2]^c), the coupling of transverse and longitudinal motions results in an unstable region of 
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Figure 3: Driving rotational transitions in a static electric field. Molecules enter the field in the 
{N, M) = (0, 0) state, and exit in the {N, 0) state. The power needed for a 7r-pulse is plotted 
as a function of the static electric field strength, for the experimental parameters given in the 
text. Curves are labelled by thek value of N. Note the logarithmic scale. 



phase-space within the separatrix, as investigated in fZTl. The accepted molecules come from a 
region close to the origin, along with the 'halo' region around this. In the transverse direction, 
part (d), molecules coming from regions close to the rods of the decelerator tend to be focussed 
too strongly at low speed and so are lost. 

In our supersonic beam, the rotational temperature is typically about 4 K, and the popula- 
tion in the (5, 0) state is only a few percent of the ground state population. To increase the 
flux of slow molecules, we will need to transfer population into our chosen state. This could 
be done by driving the microwave transition from the ground state in the presence of a static 
electric field. Without the static field it is impossible to drive this transition, since the total 
angular momentum quantum number has to change by 5 units. In the presence of a field, states 
of the same M but different N are mixed and the selection rule on N is no longer a good one. 
In the limit of very strong fields the molecules are strongly polarized along the field axis and 
the eignestates, known as pendular states, are those of a two-dimensional angular oscillator 
ir29 |. They are labelled by M and by a vibrational quantum number Vp, which only changes 
by one unit in an electric dipole transition. The state that correlates to the field-free (5, 0) state 
is the one with {vp, M) ~ (10, 0), while the ground state has (up, M) — (0, 0). It follows that 
the required transition is forbidden in the strong-field limit as well as in the weak-field limit. 
At intermediate fields, it may be possible to drive the transition. 

We suppose that ground state YbF molecules, travelling at 340 m/s, cross a resonant mi- 
crowave beam at right angles. The microwaves have a cylindrical focus, 20 mm long in the 
direction of the molecular beam, and 4 mm in the perpendicular direction. Figure [3] shows 
the power required to drive a 7r-pulse to the state that correlates to the field-free {N, M — 0) 
state, as a function of the applied, uniform, static electric field strengtlj^ As anticipated by 
the considerations above, the required power is high at both low and high fields. The power 
requirement increases rapidly with N, and the electric field that minimizes this power also in- 
creases with N . For the (5, 0) state, the power is minimized in a static field of 49kV/cm, and 
is then 2 W at 286 GHz. This seems unfeasible. Furthermore, the power requirement increases 

tXo apply this plot to a different rigid rotor molecule, X, multiply the power by (^YbF / M-x)^ and multiply the 
electric field by (Bx/-BYbF)(MYbF/Mx) 
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if the transition is inhomogenously broadened, and it will be unless the static field is uniform 
to 1 part in 10^. A better route takes the molecule via the (3, 0) state. At 31 kV/cm the transi- 
tions correlating to (0, 0) — > (3, 0) and (3, 0) — >■ (5, 0) have the same frequency, 126 GHz. At 
this field strength, the power requirements are 232 /iW and 7 /iW respectively. Inhomogeneous 
broadening is severe for the first transition whose Stark shift is about 2GHz per kV/cm at this 
field. If the field in the interaction region is non-uniform at the 10~^ level, the power needed 
increases to about 200 mW, a demanding, but not unfeasible, requirement. 

2.2 Alternating gradient deceleration in the ground state 

Whereas weak-field seeking molecules are naturally focussed through a Stark decelerator, 
strong-field seekers are naturally defocussed. The alternating gradient (AG) decelerator pro- 
vides a solution to this focussing problem, as explained in detail in 121] and ll32l . The 
molecules move through an alternating sequence of focussing and defocussing lenses, always 
passing closer to the axis in the defocussing lenses than in the focussing lenses. Because the 
transverse force is proportional to the off-axis displacement, these molecules are on stable tra- 
jectories. In the ideal case, each lens focusses the molecules in one transverse direction, and 
defocusses them in the other, the force components being — kx and Fy = —ky. The two 
directions alternate in successive lenses. The stability of the trajectories, and the transverse 
phase space acceptance, are determined by the product fir, where fl ~ \J~kJm is the angu- 
lar oscillation frequency inside a focussing lens, and r is the time spent in each lens. When 
J7r ^ 1, the focussing by one lens is almost exactly counteracted in the next and the net effect 
is very weak focussing. The trajectories are stable but have large amplitude, so the accep- 
tance is small. As fir increases, the transverse extent of the beam modulates, always being 
larger in the focussing lenses than in the defocussing lenses. Molecules with greater transverse 
speeds can now be transmitted, and so the acceptance grows. At some critical value of fir, the 
focussing becomes so strong that the molecules are overfocussed and all trajectories become 
unstable. 

AG deceleration of ground state YbF has been demonstrated experimentally li20J using 
a short prototype machine. Other strong-field seeking molecules have also been decelerated 
using similar machines lfT9l [30l [3T1 . In these experiments, the degree of deceleration has been 
limited and the transmission rather low. The AG decelerator is particularly sensitive to lens 
aberrations which tend to reduce its transmission enormously. Here, we consider a particular 
design where we aim to minimize the detrimental effects of aberrations and provide a very 
large reduction in the velocity. The design is illustrated in Fig. [4]; a). Each lens is formed by a 
pair of rods, radius R, their axes parallel and separated by 2(i? + tq), similar to the geometry 
used in all AG experiments to date. Specifically, we take tq = 2 mm and i? = 9.3 mm. At the 
exit end, each rod is terminated by a semi-ellipsoid of radii 9.3 mm along x and y, and 13.3 mm 
along z. At the entrance, the edge of the rod is rounded with a 2.7 mm radius of curvature. The 
total length of each electrode is 29.3 mm, and the distance from the entrance of one lens to that 
of the next is 34.7 mm. 

The choice of i?/ro is determined by a trade-off between the strength of focussing and 
the size of the non-linear contributions to the force |f33l. When i?/ro is large, the focussing 
will be weak and the acceptance will be small. When it is small the beneficial effect of the en- 
hanced focussing is outweighed by the detrimental effect of large nonlinear forces. Our choice, 
i?/ro = 14/3, is close to optimum for a decelerator where the length of the lenses is approxi- 
mately equal to the length of the gaps between lenses, and where fir has its best value of about 
0.7. As shown in Fig. 11 of lf33l . the transverse acceptance is then 0.012 x 2/iofii?o''o/™- For 
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Figure 4: Alternating gradient deceleration of ground state YbF. (a) Projections of the decel- 
erator onto the yz- and zy-planes. (b) The z-dependence of the Stark shift, with the position 
of the central electrode indicated in grey. The entrance end of this electrode is positioned at 
z = 0. (c) Longitudinal phase-space acceptance calculated from the potential and the turn-on 
and -off positions shown in (b). (d) The focussing (dashed line) and defocussing (solid line) 
force constants, — —dF^/dx and ky = —dFy/dy. For ease of comparison, —ky is plotted 
along with +kx- 



our parameters, this is 70 mm^ (m/s)^, a factor of 2 smaller than the transverse acceptance es- 
timated for the WF decelerator above. A higher transverse acceptance can be achieved using a 
more complicated electrode structure to form each lens, but only by a factor of about 2.5 ll33l . 

The deceleration regions are the ellipsoidal ends of the rods. Figure ^h) shows the z- 
dependence of the Stark shift for our lens structure, along with a cross-section of an electrode 
so that its position relative to the map is clear. The lens is charged to ±40 kV, and its neighbours 
are grounded. To make a first estimate of the longitudinal acceptance, we consider only those 
molecules travelling along the axis of the machine, ignoring any coupling between transverse 
and longitudinal motions. Then, the position z of a molecule relative to that of the synchronous 
molecule is given by the mean (relative) acceleration ll32l 

^ Wjz^n + S)- Wjzon) - Wjz.B + S) + Wjzos) 
dt^ mD ' ^ ' 

where W{z) is the z-dependent Stark shift shown in Fig.Qb), Zon and Zos are the positions 
of the synchronous molecule when the fields of the lens are turned on and off, and D is the 
lens-to-lens distance. By solving this equation of motion using the turn-on and turn-off posi- 
tions indicated, we obtain the phase-space trajectories of non-synchronous molecules shown in 
Fig.Qc). The outer trajectory is the separatrix between stable and unstable motion. Using this 
model, we obtain a longitudinal acceptance of 270mmm/s, almost 12 times as large as that for 
deceleration in the (5,0) state (Fig.|2]). In the Az direction, the accepted region is about four 
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Figure 5: Phase-space acceptance of the first section of AG decelerator described in the text, 
projected onto the (a) (x, v^), (b) {y, Vy) and (c) {z, Vz) planes. The section has 56 stages and 
reduces the speed from 340 m/s to 266 m/s. 



times larger and is a direct result of the larger spacing between deceleration stages. In the Av^ 
direction, the accepted region is about three times larger, because of the larger Stark shift in 
the ground state. Within the model of uncoupled longitudinal and transverse motion, the total 
6D acceptance of the AG decelerator is 18, 900 mm'^(m/s)'^. 

In reality, the transverse and longitudinal motions are rather intimately coupled. We will 
arrange the lenses so that fir is close to the optimum value for the synchronous molecule. 
The non-synchronous molecules will have different values of this parameter. Referring to 
Fig.Qb), a molecule that is ahead will spend less time in the lens and so will not be fo- 
cussed as effectively, while a molecule that is behind will spend more time in the lens and 
may be overfocussed and become unstable. Furthermore, deceleration is always accompanied 
by transverse defocussing. This effect is shown in Fig.Qd), where the transverse force con- 
stants kx = —dFx/dx, ky = —dFy/dy are plotted as a function of z. Inside the lens, where 
there is no longitudinal electric field gradient, these two force constants are equal and opposite, 
kj. = —ky = 11 aN/m. Moving into the fringe fields of the lens, the focussing force constant 
decreases along with the electric field, whereas the defocussing force constant increases to a 
value that depends on the longitudinal curvature of the electrodes. The smaller the radius of 
curvature, the larger this defocussing peak becomes. We have chosen a fairly gradual termina- 
tion of the electrodes at the exit end, where the deceleration occurs, and the magnitude of ky 
peaks at 15 aN/m. At the entrance end, the electrodes begin suddenly, and so the defocussing 
constant rises to a much larger value of 25 aN/m. However, molecules that stay close to the 
synchronous phase never see this excess defocussing, because the voltages are not yet on when 
they pass through this region. Since the focussing is most effective when the gaps between suc- 
cessive lenses are kept short, it is good to keep this entrance region short; that is why we have 
opted for an asymmetry between the entrance and exit ends. When the turn-on and turn-off 
positions are as indicated in part (b) of Fig.|4] the mean value of ky seen by the synchronous 
molecule exceeds that of k^ by 19%. This makes the trajectories unstable if fir is too small, 
making it even more important to keep this parameter close to its optimum value. Provided 
that is done, the excess defocussing has a rather small effect. We note that molecules lagging 
behind the synchronous one see the very large defocussing constant at the lens entrance, and 
we can expect these to be lost from the decelerator. 

If the turn-on and turn-off positions are kept constant, the value of r increases as the 
molecules slow down. Eventually, the limit of stability is reached and the molecules are lost. 
To avoid that, the lens length seen by the molecules needs to be reduced in harmony with the 
deceleration: the lenses should be long at the beginning and short at the end. Since deceleration 
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occurs at the end of every lens, it is better to use n short focussing lenses, followed by n short 
defocussing lenses, instead of using long lenses and alternating each time. By using a larger 
value of n at the beginning than at the end, efficient deceleration can be combined with stabil- 
ity of focussing. We have simulated the motion of ground state YbF molecules through an AG 
decelerator having the electrode structure described above. The decelerator was divided into 
four sections, with n = 4, 3, 2, 1 and with 56, 36, 32, 15 deceleration stages in these sections 
respectively. The first section reduces the speed from 340 m/s to 266 m/s, and the subsequent 
sections reduce this further to 205, 128 and 66 m/s respectively. In order to avoid a mismatch 
between the phase-space distribution exiting one section and the acceptance of the next, the 
turn-on position was gradually adjusted within each section so that the value of fir remained 
approximately the same throughout the deceleration process. 

Figure |5] shows the accepted phase-space volume of the first section, projected onto the 
{x,Vx), {y,Vy) and {z,Vz) planes. The plots were obtained as described above for the WF 
decelerator. In the {x,Vx) plane the accepted molecules form a diverging beam because the first 
four lenses of the decelerator focus in this direction, while in (y, Vy) they form a converging 
beam because those lenses defocus in that direction. In the {z,Vz) plane, the shape of the 
accepted area is similar to the calculated one shown in Fig.^c), though a little smaller In 
each plane, there is a dense central core of accepted molecules, with a less dense region near 
the edges. These regions of lower density are due to the coupling between the two transverse 
motions, and between the transverse and longitudinal motions. The acceptance of this first 
56-stage section is 2750 mm^(m/s)'^, more than twice that found for the first 200 stages of the 
WF decelerator, where the final speed was nearly the same. Taken alone, the acceptances of 
the other three sections are similar, but when we concatenate them, the net acceptance falls. 
We calculate 850 mm'^(m/s)'^ for the first two sections, 340 mm'^(m/s)^ for three sections 
and 254mm^(m/s)'^ for all four There are several causes for this additional loss. Firstly, 
those molecules near the edges of the acceptance region in Fig.|5]are on metastable trajectories 
and will be lost in subsequent sections. Secondly, there remains some mismatch between one 
section and the next, despite our efforts to control this. Thirdly, as the turn-on position is 
moved towards the turn-off position the excess defocussing at the end of each lens plays an 
increasingly important role and can result in additional loss. Nevertheless, the acceptance of 
the AG design studied here is similar to that of the WF decelerator and it is likely that further 
design iterations will improve on this. 

This 139-stage AG decelerator removes 96% of the initial kinetic energy, but the molecules 
are still not slow enough to be trapped. Five more stages are needed. These final stages are 
particularly difficult to design because of the large fractional changes in speed, and they would 
need to be designed in conjunction with the trap to ensure that the molecules remain focussed 
and are coupled efficiently into the trapping region. An alternative is to use a combination of 
AG and WF deceleration. In this scheme, the AG decelerator would be used to remove the 
majority of the energy, and then the molecules would be switched into the weak-field seeking 
state for the final reduction of the velocity. 

3 Trapping 

The electric field of the trap is also the field that interacts with the edm being measured. The 
edm signal is proportional to the integral over time of the polarization factor, P, which in turn 
depends on the electric field as shown in Fig.[T|c). The trap should have a bias field so that 
\P\ remains large and approximately constant, i.e. close to its turning point for a weak-field 
seeking state, or close to unity for a strong-field seeking state. 
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It is not difficult to design a very deep, suitably-biased, electrostatic trap for weak-field 
seeking YbF molecules in the (5, 0) state. As an example, we consider using the "chain-link" 
trap discussed in |34|. This has a minimum field at its centre, which we choose to be 130 kV/cm 
in order to maximize the value of |P| and minimize its variation. For modest applied voltages, 
the phase space acceptance of such a biased Stark trap is of order 10'' mm^(m/s)'^, very much 
larger than the phase-space volume occupied by the molecules coming from the Stark deceler- 
ator So, with a suitably careful coupling of the molecules into the trap, all the available slow 
molecules could be trapped and used in an edm measurement. 

If we are instead to trap molecules in the strong-field seeking ground state, an ac trap is 
required. For this, we consider a cylindrical ac trap of the kind first proposed by Peik |35 1 and 
recently used to trap ammonia molecules |36, 37 1. The centre of the trap is a saddle point, 
and the confining and deconfining directions alternate so that the molecules are trapped dy- 
namically - the same principle that is used to guide molecules through the AG decelerator. An 
electric field of lOOkV/cm at the centre of the trap would gives a polarization factor of 0.86. 
For the ammonia experiments, the acceptance of this trap was found to be 270 mm'^(m/s)^. 
Scaling to the YbF case, with its larger mass and larger dipole moment, we obtain an ac- 
ceptance of about 50 mm^(m/s)'^. This small acceptance makes trapping in the ground state 
rather unattractive for an edm measurement. Additionally, the large currents associated with 
the switching of the high voltages are likely to be incompatible with the extremely high level 
of magnetic field control needed for the experiment. 



4 Statistical sensitivity 

We now have all the information needed to estimate the statistical sensitivity of a trap experi- 
ment relative to that of our beam experiment, using Eq. ([TJ. It is convenient to write the total 
number of participating molecules as = pVrT, where p is the phase-space density, V is 
the phase-space acceptance of the experiment, r is the mean number of shots per second and 
T is the total integration time. The first and last factors will be common to both beam and trap 
experiments and so factor out in the comparison. We introduce a figure of merit which is the 
ratio of the sensitivity in a trap experiment to that in our current beam experiment, 

g ^'^ejbcam |-^trap| ^trap / ^rap / ^trap ^<y^ 

^'^e.trap |-^bcain| 'Hocam V ^^cam y ^bcain 

In the beam, the volume of phase space occupied by the participating molecules is Vbcam ~ 
6 X lO'* mm'^(m/s)^. The other parameters are rbcam = 12.5 s^^, |i-bcam| =0.7 and Tboam = 
10-3 s. 

Consider first an experiment where the number of available molecules is limited by the ac- 
ceptance of the WF Stark decelerator modelled above, which we found to be 350 mm'^(m/s)3. 
We take a coherence time in the trap of 1 s, and we suppose that the trap is reloaded every 
1.2 s. With a polarization factor of 0.36 we obtain a figure of merit of S ^ 10. This would be 
a very significant improvement in the sensitivity of the experiment, offering the possibility of 
reaching well below lO^^^e.cm in statistical uncertainty. Note that the repetition rate of the 
proposed trap experiment is limited by the coherence time, not the pulse rate of the source, and 
so its sensitivity could be improved even further by loading multiple traps, though at a cost of 
further experimental complexity. 

Finally, let us consider an experiment that is limited by the acceptance of the trap it- 
self rather than the method used to fill the trap. This may become possible by loading the 
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trap from a buffer gas source, or from an improved decelerator design. In this case Vtiap ~ 
lO"* mm'^(m/s)'^ and we find S « 50. All other factors being equal, in particular the phase 
space density produced at the source, the statistical sensitivity would then be about 2 x 10~^^ e cm 
after T — 24 hours of integration. To realize such exceptional sensitivity we would, of course, 
need to reduce all other sources of noise to this level. We see no fundamental reason why that 
should not be possible. 

5 Physics in the trap 

So far we have considered how the statistical sensitivity of an edm experiment might be im- 
proved by decelerating YbF molecules and loading them into a trap. Now we consider the 
impact of the trap environment on such an experiment. Our present edm experiment measures 
the precession of the molecular spin in a combination of nominally parallel electric and mag- 
netic fields. We prepare the spin in a direction perpendicular to the fields and, after a fixed 
time, we measure the accumulated precession angle, which depends on the magnetic moment 
and the electron edm. In the trap, the electric field is necessarily inhomogeneous - that is how it 
exerts a force. Therefore, as molecules move around within the trap, they experience a substan- 
tial variation in the magnitude and direction of the field and this fluctuating environment varies 
from one molecule to another. It is therefore essential to consider what happens to the spin 
precession angle when the direction of the local electric field at the position of the molecule 
is changing with time. In the following, we consider two consequences of this: an effective 
magnetic field inhomogeneity, and a geometric phase. 

5.1 Effective magnetic field inhomogeneity 

Because of the large tensor Stark shift induced by the electric field of the experiment, only 
the magnetic field component parallel to the local electric field direction contributes to the 
precession angle 1, 12j . The accumulated phase angle due to a magnetic field B is 

^B = ^ jBiv{t))-e{f{t))dt, (8) 

where g and are the g-factor and Bohr magneton, e is a unit vector along the local electric 
field direction, and the integral is taken over the trajectory of the molecule, f(<). 

In a biased electric trap (Sec.[3]l, the trapping field is superimposed on a uniform bias 
electric field. We define z = e{0), the electric field direction at the centre of the trap. For 
those molecules that stay close to the centre, the local electric field direction is always close to 
z, whereas molecules that reach the outer regions of the trap will see e making larger angles 
with z. It follows that the integral in Eq. ^ changes from one trajectory to the next even if the 
magnetic field is perfectly uniform. This spread in the accumulated phase is a source of noise 
in the experiment and may degrade the experimental sensitivity if it is too large. The size of 
the phase spread depends on the temperature of the molecules relative to the Stark shift at trap 
centre (expressed in temperature units). If the temperature is high the molecules will explore 
a large range of electric field directions and the spread will be high. If cold, the molecules 
remain close to the trap centre where the electric field direction is close to z everywhere. 

To explore this, we simulated the motion of molecules in the chain-link trap discussed in 
Sec.|3] and evaluated 0b for each trajectory using Eq. (j8]l. The molecules were drawn from an 
initial phase space volume of (4 mm x 4m/s)'^. As discussed above, the molecules are most 
sensitive to fields along z and so we applied a uniform magnetic field, B^, in this direction. We 
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find that both the mean and the standard deviation of the accumulated phases increase linearly 
with the coherence time, t. Specifically we find a mean of {4)b) ~ 0.9 0o and a standard 
deviation of A0s ~ 0.05 4)q, where 0o = giJ-sBzT /his, the phase accumulated by a molecule 
that remains at the centre of the trap throughout. This phase spread should be compared to the 
error in a shot-noise limited measurement of the edm phase, which is 1/VN rad (see Eq. 
Therefore, to ensure that the magnetic field does not increase the noise and so degrade the 
sensitivity, we require /S.(f>B < 1- For r = 1 s, the corresponding requirement on magnetic 
field is Bz < 200 pT. Magnetic field gradients will also contribute noise as different molecules 
sample different regions of the trap. Our simulations show that, for a coherence time of 1 s, 
these gradients need to be smaller than 50 nT/m if they are not to degrade the sensitivity. These 
field requirements are not too demanding; they can be satisfied using a multi-layer magnetic 
shield with an overall shielding factor of 10^, e.g. Il38l . 

5.2 Geometric phase 

We show below that, for a sufficiently strong electric field, the motion of the spin follows 
the electric field direction adiabatically, and that the spin polarization will always lie in the 
plane perpendicular to the local electric field. In this adiabatic limit, one might guess that 
the accumulated phase of the spin would be the time-integral of the instantaneous precession 
frequency. However, in the case of a rotating magnetic field, it is known from the work of Berry 
ll39l that the accumulated phase angle has a second component, known as the geometric phase. 
This depends on the history of the field direction but not on the precession dynamics. We might 
expect that a similar geometric phase will appear in the spin precession of a molecule when 
the electric field changes direction. If we are to obtain the electron edm from the precession 
angle, we must be able to control any such geometric phase with high precision. 

While the geometric phase of a spin precessing in a rotating magnetic field has been stud- 
ied in great detail 1391 . the case of a rotating electric field is less well-studied and the few 
treatments that do exist tend to be phrased in the somewhat arcane language of differential ge- 
ometry [40J. In the interest of clarity, therefore, we first present a simple quantum mechanical 
treatment of the geometric phase in a rotating electric field, before going on to consider the 
impact of this phase on a potential edm experiment. 

In the limit of large electric field, which is well-fulfilled in our experiments, the magnetic 
field makes no significant contribution to the geometric phase evolution. This is because the 
large tensor Stark shift strongly suppresses the interaction with magnetic field components 
perpendicular to the electric field 11 2J . Therefore we neglect the applied magnetic field in the 
following treatment, even though it is possible in principle to include it BTIl . 

5.2.1 Spin-one system in a rotating electric field 

Let us consider a spin-one system, corresponding to the F — 1 state of the YbF molecule used 
in our edm experiment. The Stark interaction mixes many states of the molecule, but for our 
present purpose it is adequately described by an effective Hamiltonian restricted to the F = 1 
manifold, consisting of a scalar and a rank-2 tensor With the electric field directed along the 
z-axis this takes the form 

HE^=S{E)f^"^ +e{E)f^^\ (9) 

where the coefficients S{E) and e{E) are phenomenological parameters that depend on the 
electric field strength E. If desired, they can be derived from a full calculation of the Stark 

" (n) 

shift including all of the molecular levels. Tg ' is the z-component of the rank-n irreducible 
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spherical tensor operator. The matrix elements of this operator in the \F = l^mp) basis can 
be written, in order of ascending m^, as 



^0 0\ 

He.^\0 a , (10) 
^0 0/ 

where we have simplified the parametrization by redefining the zero of energy and introducing 
A = ^JY/2,5{E) — ■\/2/15 e{E). We have dropped the explicit functional dependence of the 
parameter A on the electric field in order to streamline the notation. 

When the electric field rotates to a new direction given by the Euler angles {a(t) , , 7(t) }, 
the Hamiltonian becomes 

H{t)^R{t)-^HE^R{t), (11) 

where R{t) is the rotation operator We adopt the Euler angle convention of Weissbluth [42 1: a 
rotation about the z-axis by 7, a rotation about the y'-axis by /3, and finally a rotation about the 
2;"-axis by a. The rotation matrix acting on states and operators (rather than on coordinates) is 
then given by 



R{t) = 



^ e*("+T)cos2 (f) ^e*"sin/3 e*("-T) sin^ (f ) ^ 

— ^e'''sin/3 cos/3 ^e~*'''sin/3 
^ eK''-")sin2 -^e-*" sin/3 e-*("+'^) cos^ ^ y 



(12) 



We note that the cylindrical symmetry of He- around z also makes H{t) symmetric around 
z" and therefore independent of aj^jThe Schrodinger equation takes the form 

R{t)-^HE^R{t)C{t) = ^h—C{t) , (13) 

where C{t) is the column matrix of coefficients Cmp from the expansion \ip) = Cmp (t) |to_f). 

The same state can be expanded on the basis that rotates with the electric field and has z" 
as its quantization axis. This has expansion coefficients Cii{t) = R{t)C{t). In terms of these. 



Eq. (13 I becomes 



R{t)-'HEM{t)R{t)-'CR{t) - th^^R{t)-'CR{t) , (14) 



which simplifies to 



where 



[He^ +hG{t)) Cnit) = ih^CR{t) , (15) 
G{t)^~iRit)-RityK (16) 



tExplicitly, 



^ I v^Asin^^ -e-'TAsin2/3 - v/2e-2»T A sin^ ^ 

H{t) = I -e*TAsin2^ 2-/2Acos2/3 e-'TAsin2;3 

\ -V2e2'TAsin2/3 e'TAsin2/3 \/2Asm2/3 
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We see that the evolution of the state expressed in the rotating z" basis is governed by an effec- 
tive Hamiltonian made up of the static II e^, plus a term hG(i) involving the time-dependence 
of the rotation. This second term is much like a fictitious force in classical mechanics. As we 
shall see, it generates a geometric phase. 



We can write G(t) explicitly, using ( 16 1 and ( 12 1: 



-\/2(d + 7COS/3) e*"(7sin;9 + i/3) 
G(t)^^ \ e-*"(7sin/3-i/3) e'" (7 sin /3 + | , (17) 

e-*"(7sin/3- y2(d + 7Cos/3) 

where dotted quantities are derivatives with respect to time. We wish to consider the adiabatic 
limit, where the tensor Stark splitting A is always much larger than the rotation rate. In this 
limit, it is a good approximation to replace G(t) by 



Gait) = . (18) 




The subscript a indicates that this is the effective operator in the adiabatic limit. This approxi- 
mation follows from the formal properties of matrices that the eigenvalues and eigenvectors of 

, a h Q . ,aOO, 

&* A f) j approach those of the diagonal matrix ( A j in the limit of large A. It is also 

familiar from time-independent perturbation theory, which applies in the limit of slow rota- 
tions, where the effects of the off-diagonal elements are smaller than those of the diagonal 
elements Ga by a factor of order Gij/ IS.. 

The operator Ga [t) is the central result of this section: in the adiabatic limit, the dynamics 
of a spin-one system in a rotating electric field are the same as those in the static field plus 
a time-dependent fictitious magnetic field directed along the electric field. This breaking of 
degeneracy between the mp = ±1 sub-levels has nothing to do with the static interaction. On 
the contrary, it is a vector property caused by the vector nature of rotatior{§] We also note that 
the adiabatic condition only requires the splitting A between the mj? = and toj? = ±1 sub- 
levels to be large compared with the rotation rates. The splitting between mp = ±1 sub-levels 
is unimportant as the angular momentu m op erator does not couple these states. 



In order to calculate Ga from Eq. (18 1, we need the Euler angles as a function of time. 
The electric field rotation fixes /3{t) and ■y{t) unambiguously, but we have complete freedom 
to choose a{t) because of the cylindrical symmetry of the Stark interaction around z". It is 
convenient to pick a{t) such that Ga vanishes, which requires 

a(i) =ao- / 7(T)cos(/3(r))dr , (19) 
Jo 

where ao is a constant of integration. In this particular rotated basis, the adiabatic Hamiltonian 
is just that of the molecule in a static electric field. However, when the field makes some 



excursion before returning to its original direction, the final value of a, given by ( 19 1, will not 
be the same as the starting value ao- This is the essence of the geometric phase. 

Figure[6]illustrates the geometry of our problem. A unit vector along z" follows the rotation 
of the electric field, mapping out the path F on the surface of a sphere. At any given point on 
the path, with coordinates 7), the x" axis lies in the tangent plane, and makes an angle a 
to the great circle of constant 7. Following an initial choice of ao, the subsequent evolution 



^Rotation is generated by angular momentum, a vector quantity. 
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Figure 6: Geometry of the problem. The z" axis is along the electric field direction defined 
by (/3,7). When the field rotates, a unit vector along z" follows the path F (red line) on the 
surface of a unit sphere. The x" axis is tangent to the sphere and lies at an angle a to the 
dashed great circle. This angle changes under parallel transport around the loop, as given by 



Eq. ( 19 1, and does not in general return to its original value oiq at the end of the loop. The 



change in a is equal to the solid angle subtended at the origin by the red curve F. 



of a under Eq. ( 19 1 is given by parallel-transport of the x" -y" plane along the path F. This 
means that for each infinitesimal displacement along the path, the new basis is the one most 
nearly parallel to the old basis. In the interest of brevity, we simply state this here and do not 
prove it ||431 . At the end of a rotation around a closed path F, it is easy to verify that a — a^, 
as given by Eq. ( 19 1, is equal to the solid angle Vl subtended at the origin by the curv^ This 



is an example of a more general result from differential geometry, that the so-called holonomy 
angle introduced by parallel-transport around a curve is equal to the integral of the Gaussian 
curvature over the region bounded by the curve. For a sphere, which has constant Gaussian 
curvature l/r^, this integral is equal to the solid angle. 

This phase shift has a mechanical consequence. Consider a spin prepared at i = along 
X = x"(0). Under adiabatic rotation of the electric field, the spin undergoes parallel transport 
into the new x" direction. At the end of a closed loop, although the spin has returned to the 
initial x — y plane, its direction x"{t) is now rotated through an angle a(t) — ao = 51, equal 
to the solid angle subtended by the path F. It is this spin rotation that poses the main challenge 
for measuring the electron edm using trapped YbF. 



^More correctly, it is f2 — 27r, but the 27r comes from our definition of the Euler angles and has no physical 
consequence 
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Figure 7: x-y projection of the normalised local electric field at the positions of 32 molecules 
evolving in the trap for 10 ms. Each circle shows the curve swept out on the sphere by the tip 
of the electric field direction vector for a single molecule, as viewed from the north pole. 

5.2.2 Geometric phase in the trap 

We have made a numerical simulation of thirty-two molecular trajectories in the chain-link 
trap of Sec.[3]in order to evaluate the effect of the geometric phase. This calculation was made 
as realistic as possible by using finite-element-analysis to determine the field. The molecules 
were loaded into the trap with a uniform distribution over ±2 mm in all three spatial directions 
and over ±2m/s in all three velocity components. Figure [7] shows x-y projections of the unit 
vector along the electric field at the position of each molecule. These views, seen along the 
the z-axis, show the path T (illustrated schematically in Fig. (|6]l) that determines the geomet- 
ric phase. It is immediately clear without further calculation that the accumulated geometric 
phases are large and different for each molecule. The solid angle for each orbit of the trap 
can approach a sizeable fraction of 27r. These simulations last for 10 ms, during which time 
the molecules execute 3 or 4 oscillations. The result is a large inhomogeneous spread of the 
spin direction in the x-y plane that completely depolarizes the molecular ensemble. Therefore 
the spin decoherence time is only a few ms, not the 1 s proposed in Sec. |4] The geometric 
phase is clearly the primary obstacle to making an edm measurement, or indeed any other spin 
precession measurement, in an electrostatic trap. 

6 Conclusions and outlook 

The use of trapped polar molecules to measure the electron edm offers a very substantial in- 
crease in statistical sensitivity. By loading an electrostatic trap using a Stark decelerator a 
sensitivity gain of 10 is feasible. Deceleration is possible in either the ground state or a ro- 
tationally excited state with comparable efficiencies. The efficiency of AG deceleration is 
particularly sensitive to the electrode geometry and timing sequence and it seems likely that its 
losses could be reduced further through refinement of these parameters. A buffer gas source 
might be able to deliver more molecules to the trap and so increase the statistical sensitivity 
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even further. 

However, the trap is a difficult environment for an edm experiment. The spatial variation in 
the direction of the electric field which provides the trapping force has two deleterious effects. 
The first is an inhomogeneous broadening of the magnetically-induced spin-precession angle, 
which can be kept under control by making the applied magnetic field small enough. The 
second is an inhomogeneous geometric phase. The long coherence time that seemed to be 
possible in the trap is precluded by this geometric phase which scrambles the spin precession 
signal in just a few milliseconds. 

One possibility for overcoming this problem might be to use a molecule with fi-doublet 
structure, for example ThO |44|. In these molecules there are states with equal mp but op- 
posite sensitivity to the edm. A superposition of two such states prepared within the same 
molecule would experience a common geometric phase shift but a differential edm phase shift. 
This might provide the basis for a reduced sensitivity to the geometric phase, but it would be 
vulnerable to differential Stark shifts between the two il-doublet levels resulting from the vari- 
ation in electric field strength. Further analysis is required to determine if this could be used 
as the basis for a workable edm measurement. A molecule with a pure spin 1/2 ground state, 
and therefore no tensor Stark splitting, would not adiabatically follow the electric field and so 
would be immune to the geometric phase arising from the rotating electric field. However, 
the molecule would lose its anisotropic response to magnetic fields which is a key feature for 
suppressing both the magnetic geometric phase and the troublesome v x E effect |8 1. Another 
possibility is to use the decelerator as a source for a molecular fountain where it would seem 
to be much easier to control the electric and magnetic fields and thereby control the systematic 
effects. The coherence time approaches 1 s in a 1 m-high fountain. Note that the molecular 
beam exiting the decelerator would need to be expanded to «10cm diameter to obtain the 
required degree of collimation. 

The difficulties associated with the trap environment occur because the temperature of the 
molecules is similar to the trap depth and so they explore most of the trap. If the molecules 
were far colder they would remain near the trap centre where the field is very uniform. This 
might be achieved by loading from a source with much higher phase-space density, so that 
we could afford to discard all but the coldest fraction, or by actively cooling the molecules 
to ultralow temperatures. There has been rapid progress in the formation of deeply -bound ul- 
tracold polar molecules using photoassociation, Feschbach resonance, and coherent transfer 
techniques Il45ll46ll47l . and these methods might be extended to molecules suitable for edm 
measurement. Several research groups are now exploring the direct application of laser cooling 
techniques to molecules 1481 . In this context we note that the vibrational Franck-Condon struc- 
ture of YbF P9l is favourable for cycling transitions. Other routes to ultracold polar molecules 
include sympathetic cooling (SU] and cavity-assisted cooling 151]. If a scheme can be devised 
to cool a suitable heavy polar molecule to low temperature, edm measurements would not just 
be improved but revolutionized. 
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